Relative performance of ancilla verification and decoding in the [[7,1,3]] Steane code 
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While ancilla post-selection is a common means of achieving fault-tolerance in quantum error- 
correction, it can lead to additional data errors due to movement or hold operations, and alternatives 
to post-selection may achieve lower overall failure rates due to avoiding such errors. We present 
numerical simulation results comparing the logical error rates for the fault-tolerant [[7,1,3]] Steane 
code using the technique of ancilla verification vs. the newer method of ancilla decoding, as described 
in [J]. We simulate a realistic QEC procedure in which failed ancilla creation requires storing the 
data until a new ancilla can be created; we find that the decoding method, which avoids the need 
for such storage, has an advantage when the failure probability is sufficiently high. We additionally 
analyze the effect of different classes of physical error (initialization, measurement, hold etc.) on 
the relative performance of these two methods. 



INTRODUCTION 



Ancilla postselection in quantum error 
correction 



Quantum error-correcting codes (QECCs) provide a 
means to protect quantum data against noise, by en- 
coding quantum states into larger Hilbert spaces such 
that some class of error operations are correctable @, H| ■ 
In a physically-realistic system, however, one must take 
into account that the correction operations themselves 
must be implemented using imperfect quantum opera- 
tions. This leads to the common requirement that for 
quantum error correction (QEC) operations to be effec- 
tive in protecting the data, they (and other operations on 
the data) must be implemented in a fault-tolerant way 
Q i.e. in a way such that a single faulty quantum gate 
cannot lead to multiple errors on the data. 

Some logical operations are naturally fault-tolerant, 
such as transversal gates in which every physical qubit 
is acted on by a separate gate, and the absence of any 
" 'cross-talk" ' between different physical qubits in an en- 
coded block means there is no opportunity for errors to 
spread. In other cases, however, ancillary states which 
interact with the data as part of the QEC process need 
to be created using non-fault-tolerant circuits, and fault- 
tolerance is instead enforced by post-selection (ancilla 
verification), in which only those ancillas which satisfy 
some measurement outcome after being created are sub- 
sequently used to interact with the data[f| Q. Thus it 
is not known beforehand whether a given ancilla will be 
used and one may need to perform multiple ancilla cre- 
ations before obtaining one which passes post-selection. 

This raises the question of how to ensure that a post- 
selected ancilla is available for QEC with sufficiently high 
probability so as not to significantly increase the overall 
failure rate of the QEC. Two obvious approaches are to 
either create ancillas sequentially until post-selection is 
passed, or to create many ancillas in parallel so there is 
a high probability that at least one will pass 0- Both 



approaches have somewhat analogous disadvantages: if 
an ancilla fails, sequential creation requires the data to 
wait (and, in general, accumulate errors) while a new an- 
cilla is created. Parallel creation avoids this, but has the 
problem that in a given physical architecture there may 
be a very limited number of ancillas that can be created 
close to the data. Thus, a passed ancilla may be created 
some distance away and need to be moved into contact 
with the data, both data and ancilla accumulating errors 
while this occurs. 

In this paper we consider addressing this problem using 
ancilla decoding [l[ , a technique originally devised to ad- 
dress the different problem of slow qubit measurements. 
As we discuss below, ancilla decoding removes the need 
for post-selection and thus guarantees that any created 
ancilla can be used with the data. Hence this ancilla can 
be created next to the data and QEC performed without 
any additional movement or waiting. Depending on the 
circuit layout and gate errors and timings, this may re- 
sult in lower overall logical failure rates and hence in the 
case of a compuation using a concatenated QECC, fewer 
resources required to achieve a given logical error rate. 



B. The Steane code and ancilla creation 

In this work we compare the logical error rate Pl ob- 
tained when performing fault-tolerant QEC operations 
for the well-known [[7,1,3]] Steane code 0], performing 
either ancilla verification or decoding. The Steane code 
encodes a single logical qubit into the state of seven phys- 
ical qubits, and can correct any error on a single qubit. It 
is a CSS stabilizer code, whose list of stabilizer generators 
is given in Table [TT] in the appendix. In our analysis we 
consider the Steane ancilla technique, in which the code 
stabilizers are measured by copying error information to 
ancillas in encoded logical states. 

A (non-fault-tolerant) picture of part of the QEC pro- 
cess is shown in Figure [T] in which the data (qubits \qd)) 
interacts with a prepared logical |0l) (where subscript 
L denotes logical states) ancilla state (qubits |o)) to de- 
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FIG. 1. A Steane code circuit with non-fault-tolerant Steane syndrome extraction, where |oi), 02), ...,|<J7) are the ancilla 
qubits. These ancilla will interact with each other via CNOT gates to produce the logical ancilla state (i.e. |0_l ) ) . Also, \qdi), 
\qd,2) ...,\qd,7) are the data qubits. The ancilla-data interactions are implemented by a set of transversal CNOT gates, where 
any Z error in the data will propagate to the ancilla. Ultimately, the ancilla will be measured in the X-basis to determine the 
error syndrome. 



termine the error syndrome for Z errors (an analogous 
interaction occurs with a logical \+l) to correct for X 
errors). We see, for example, that for fault-tolerance the 
|0l) ancilla, which interacts with the data as the source 
for a transversal CNOT gate, needs to be prepared so 
that a single gate failure will not cause the ancilla to have 
multiple Pauli X errors (likewise with the state and 
Z errors). Such errors would get transferred to the data, 
and can occur in the circuit shown in Figure [T] since the 
CNOT gates involved in preparing |0i) can propagate 
a single gate error to multiple qubits within the ancilla 
block. 

The standard approach for ensuring this using ancilla 
verification (which occurs between ancilla preparation 
and data interaction) is shown in Figure [H Passing the 
verification procedure is dependent on the outcome of the 
transversal Z measurement performed on the verifier. An 
analogous procedure occurs for Z errors and the \+l)- 

The alternative method of ancilla decoding [l| is il- 
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FIG. 2. An example of a single ancilla preparation and veri- 
fication circuit taken from the Steane code with Steane syn- 
drome extraction. Where \Q) is the encoded ancilla state, 
and, \V) is the encoded verifier state. Any X errors which 
have occurred in the creation of \Q) will propagate to the 
verifier via the transversal CNOT. 



lustrated in Figure [3] The basic principle is that after 
interacting the ancilla with the data, one applies a de- 
coding operation to the ancilla to extract both the error 
syndrome (as usual) and to determine the presence of 
any errors on the ancilla which may have been trans- 
ferred to the data (X or Z, depending on the ancilla in 
question). In addition, the ancilla interacts with a sec- 
ond ancilla block which is in a product of |0) states (i.e. 



3 



Data qubits 



Encoded ancilla 
state. 



Second ancilla 
block 



Decoding 
circuit of the 
first sub-block 



Decoding 
circuit of the 
second sub-block 



FIG. 3. The ancilla decoding procedure. After interacting the 
data with a non-postselected ancilla a decoding procedure is 
applied to determine the error syndrome and any multi-qubit 
errors which may have been transferred to the data. A second 
block allows one to distinguish between errors from encoding 
and decoding 
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FIG. 4. A schematic for the decoding circuit of the non- 
postselected |0l) ancilla for the Steane code. We are able to 
detect single- or/and multi-X qubit errors which propagate 
to the data during the ancilla-data interactions (there is an 
equivalent circuit for decoding the non-postselected an- 
cilla to detect the single- or/and multi- Z qubit errors to the 
data). After the decoding, the ancilla will be measured in the 
X— and Z— basis as shown in the Figure. Measurements in 
the X— basis allow us to detect the Z errors which were ini- 
tially in the data, while the Z— basis measurements will allow 
us to detect the X qubits errors due to the non-postselected 
ancilla. 



|0}® 7 ), which is also decoded after interacting with the 
data. This allows one to distinguish between errors of the 
first ancilla acquired during ancilla encoding (which will 
have been propagated to the data and need correcting) 
and those during decoding (which will not be transferred 
to the data or to the second block). The decoding circuit 
is illustrated in Figure |U 

Any propagated errors (single or multi-qubit) can then 
be corrected on the data. An important factor in the suc- 



cess of this technique is that, for the purposes of fault- 
tolerance, only "first-order" patterns of errors, which can 
be caused by a single gate failure (though they may still 
affect multiple ancilla qubits through error propagation, 
or through errors on two-qubit gates) need to be cor- 
rected. Furthermore, many error patterns are equivalent 
up to stabilizer operations. This allows the limited infor- 
mation from the decoder to be sufficient to correct any 
such errors. Secondly, different classes of errors (X and Z 
or vice versa) are respectively detected by the standard 
syndrome measurement and by the additional decoder 
syndrome. For example, a |0l) ancilla is used to deter- 
mine the syndrome for Z errors on the data, but is at risk 
of propagating X errors to the data, with the additional 
decoding giving information about the latter. The use 
of this method, therefore, means that any ancilla which 
is created (regardless of first-order errors) can be used 
in QEC, removing the need for verification. The encod- 
ing circuit is illustrated in Figure [5] in the appendix, each 
gate in that circuit being faulty and capable of producing 
an error which will propagate (cither as single or multi- 
qubit error) to the data via the transversal CNOT gates 
(as shown in Figure [T]). For example, if a single fault has 
been produced by the two output channels of CNOT— 7, 
or —8, or —9 (in Figure |HJ), then this error will propagate 
as a two-qubit error to the data. This is discussed in 
more detail in the appendix. 

As mentioned above, the original proposed advantage 
for ancilla decoding was to avoid long waits for ancilla 
verification in the case of slow measurements. Since, 
unless performing non-Clifford-group operations on our 
data, we can operate in the "Pauli frame" (merely clas- 
sically recording the necessary corrections and updating 
the stabilizer accordingly rather than actually applying 
them), there is no corresponding need for the data to 
wait for the outcome of the ancilla decoding operations. 
One could also partially avoid the problem of waiting for 
verification by simply beginning the verification well in 
advance, but even for fast measurements, this does not 
avoid any delays due to having to restart preparation if 
verification fails. In our analysis we consider the case of 
measurement operations no slower than any other gate, 
and show that decoding still gives an advantage for this 
reason, even if verification failure is rare. 

As discussed above, we are interested in the limitations 
imposed by practical QEC architecture, in which a ver- 
ified ancilla may not always be available. We simulate 
a case where ancilla creation occurs serially and a failed 
ancilla verification requires the data to wait while a new 
ancilla is created and compare the performance (in terms 
of the overall Pl) of a QEC procedure under these cir- 
cumstances to QEC using the decoding method. While 
we consider the specific case of the Steane code using 
the Steane ancilla technique, it is conjectured [l[ that a 
similar decoding procedure can be found for other CSS 
QECCs. 
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II. SIMULATION PROCEDURE 

To compare the two methods of ancilla interaction, we 
performed Monte Carlo simulations of the complete QEC 
procedure (so interaction of the data with two ancillas, 
one for the correction of each of X and Z errors) imple- 
mented using faulty gate operations. Our simulation soft- 
ware was QASM-P, software based on QASM by Cross 
. All gates were simulated using the following common 
stochastic error model for depolarising noise, a function 
of a single error probability p: 

1. Attempting to perform a data qubit identity I (hold 
operation), but instead performing a single qubit 
operation X, Y, or Z, each occurring with probabil- 
ity p/3. 

2. Attempting to initialize a qubit to |0) or |+), but 
instead preparing |1) and |— ), respectively, with 
probability p. 

3. Performing a Z-basis or A-basis qubit measure- 
ment, but reporting the wrong value with proba- 
bility p. 

4. Attempting to perform a CNOT gate, but instead 
performing a CNOT followed by one of the two- 
qubit operations I®X, I®Y, I®Z, X®I, X®X, 
X®Y,X®Z,Y®I,Y®X,Y®Y,Y®Z,Z®I, 
Z®X, Z®Y, or Z®Z, each with probability p/15. 

We considered a range of values of v below the threshold 
for Steane code of roughly 10~ 4 Q. All gates are as- 
sumed to have the same duration, with hold operations 
implemented by applying the identity gate the appropri- 
ate number of times. 

To determine Pj, for the QEC procedure, the data was 
first prepared, without error, in a logical cigenstate \0l}- 
The QEC procedure was then performed using the er- 
ror model above. Finally, a second QEC was performed 
with errors turned off, to project the output state into 
a logical eigenstate. This eigenstate was then compared 
against the input and a logical error reported if they did 
not match. The overall procedure is illustrated in Fig- 
ure [5] From the Monte Carlo simulation we obtain Pl 
as a function of p. Note that due to the use of the |0l) 
state, only logical X errors are detected; this suffices to 
compare the techniques since by the symmetry of the un- 
derlying physical error model we expect the rate of logical 
Z errors to be approximately the same (up to some slight 
asymmetry due to performing one correction before the 
other). 

When simulating ancilla verification, we assume a se- 
rial creation process for each of the two ancillas used in 
the QEC. We assume the ancilla is initially prepared an 
appropriate length of time prior to the QEC so that the 
data is not required to wait for ancilla preparation, so 
long as the first ancilla successfully passes verification. 
However, if ancilla creation fails, the data is held until 
a new ancilla can be created and verified. Holding the 




FIG. 5. Steps involved in the QEC simulation. 



data continues until verification is passed, at which point 
the QEC continues. We hold the data for a total of 6 
time steps (the number required to create and verify a 
new ancilla) for every failed verification, with no upper 
bound to the number of failures that can occur. Once an 
ancilla has successfully passed verification the QEC pro- 
cedure continues. To determine the extent to which the 
verification process is affected by this data holding (as 
opposed to any other element), we additionally simulate 
an unrealistically optimistic case of verification where no 
additional holding is required if verification fails. 

Additionally, we generated simulation failure rates in 
which particular classes of errors were considered in iso- 
lation i.e. with other errors turned off. As in e.g. [Ioj |. 
we define these classes of errors as follows: 

• Class 0: errors from preparation and measurement. 

• Class 1: errors from single-qubit gates (including 
the identity). 

• Class 2: errors from two-qubit gates (i.e. the CNOT 
gate). 

III. RESULTS AND ANALYSIS 

A direct comparison of the Pl for all three QEC tech- 
niques is shown in Figure |51 The results show an advan- 
tage for the decoding technique over verification over the 
whole range of errors considered, larger at larger error 
rates with a roughly 2-fold reduction as the maximum. 
Comparison with the Pl values for verification without 
additional holds, which are lower still, indicates as ex- 
pected that the increased errors in the verification tech- 
nique are due to ancilla failures, and that in the absence 
of such failures the additional operations required for an- 
cilla decoding lead to a slightly higher error rate. The 
rate of ancilla failure as a function of p is plotted in Fig- 
ure [7J Since ancilla verification detects single errors on 
the ancilla, this is roughly equal to (number of error lo- 
cations in ancilla creation circuit) xp, and the rate scales 
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FIG. 7. The ancilla failure rate vs. physical error rate for the 
[[7,1,3]] Steane code with Steane ancilla 

linearly with p as expected. Thus at low physical error 
rates ancilla failure occurs in a proportionally low frac- 
tion of QEC events, but still contributes significantly to 
the logical error rate. 

Comparing the error rates by class in Figure [8] we see 
as expected that the largest source of error in verification 
is due to hold operations, while the other two techniques 
are dominated by CNOT errors. This emphasizes the 
need to take into account nondctcrministic operations 
(such as verification) when assessing the impact of gate 
errors on QEC procedures. 

IV. CONCLUSION 

Using Monte Carlo simulations, we have compared the 
logical error rate, Pl, for implementations of the Steane 
code with Steane ancillas using two different QEC tech- 
niques: ancilla verification with serial recreation for fail- 
ure events, and the ancilla decoding procedure. We find 
that, even when measurement times are no longer than 



those for other operations and verification failures are 
rare, the decoding procedure, though otherwise more 
complex, is advantageous in avoiding data holds due 
to verification failures and leads to lower failure rates. 
Clearly these results are highly dependent on the error 
model chosen: many physical systems are dominated by 
errors in two-qubit operations, and a significantly lower 
error for holds might well result in verification being the 
preferred method. However, as well as showing that an- 
cilla decoding can be of use even in the fast measure- 
ment regime, the results demonstrate that considering 
additional detail in practical implementations of QEC 
procedures (such as the need for ancilla recreation if ver- 
ification fails) can make a significant difference to over- 
all error rates and the preferred means of implementing 
a QEC. The particular case of verification is of interest 
given that many QECCs rely on post-selection to achieve 
fault-tolerance. We note that an analysis of the effect of 
decoding in trapped-ion-based quantum computing has 
recently been published [ill]. 

One motivation for assuming a serial rather than par- 
allel ancilla creation process in the verification model was 
that we were not required to make any particular assump- 
tions about the layout of the underlying physical system. 
However, parallel creation, as a discussed in the introduc- 
tion, is an obvious possible solution to avoid verification 
delays, and for a given layout it may well be that any 
additional errors due to qubit movement in the parallel 
case are less damaging than those incurred through serial 
creation, and possibly less than those from the additional 
operations in decoding. Comparing the three techniques 
over a range of likely layouts and error models would 
therefore be a promising avenue for future research. 
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Appendix: Details of ancilla decoding 

In this section we describe the possible ancilla errors 
due to a single fault in the encoding circuit, and the cor- 
responding error syndromes when using the ancilla de- 
coding technique. This is an expansion of the summary 
given in We will consider the case of X errors due 
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FIG. 8. Logical error rate Pl vs. physical error rate p for class — errors (black circle points), class — 1 errors (red triangle 
points), and class — 2 (green square points) for (a) Verification (without additional holds), (b) Verification and (c) Decoding. 



to faults in the encoding circuit for |0l); the case of Z 
errors when encoding \+l) is directly analogous. 

Figure IH1 illustrates the circuit for encoding \0l)> which 
possesses two types of gates: (a) single qubit gates de- 
fined as preparation of |0) and |+) with waiting gates rep- 
resented by bold line segments in Figure El as well as (b) 
two qubit gates defined as CNOT gates. The decoding 
circuit can uniquely identify any single-qubit error. T.,„ 
qubit-errors caused by a single fault will either be errors 
on both outputs of a CNOT gate, or a single error which 
at some point propagates to two errors via CNOT, thus it 
suffices to further consider only the case of a CNOT gate 
on which both outputs have X errors. These are listed in 
Table Q] according to the encoding circuit labelling given 
in Figure [SI 

As shown, dual output X errors propagate to single 
X errors to the data in the cases of CNOT gates 4, or 
5, or 6; or as two X errors to the data in the cases of 
CNOT gates 7, or 8, or 9. Sec Table H However, for 
CNOT gates 1, or 2, or 3, dual X errors in the two output 
channels will propagate as four X errors on the data. 
These are equivalent to the X-stabilizer generators of the 
Steane Code, which accordingly do not affect the data. 
The stabilizer generators of the Steane Code are listed 
in |TTJ Finally, for CNOT gates 4, or 5, or 6; the single 
faults which are produced in the two output channels 
will propagate as a dual X errors on the data. 

These X errors (single and double X errors) on the 
data can be identified by the decoding circuit (Figure 0} 
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FIG. 9. Encoding circuit for \0l) 



TABLE I. Examples of output errors produced by X errors 
on both outputs of CNOT gates in the \0l) encoding circuit 
(Figure El). 



CNOT Multi-qubit data errors 


1 


V1V3V5V7 = Gi 


2 


X4X5X6X7 = G2 


3 


X2X3X6X7 = Gj, 


4 


X1X5X7 = G1X3 


5 


V4V5V7 = G^Xq 


6 


X^X-iXf, = G3V7 


7 


X2X3 = G-iX&Xj 


8 


X4X5 = G^X^Xi 


9 


X1X7 EE G1X3X5 
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TABLE II. The stabilizer generators of the Steane code. 



X-stabilizer Z-stabilizer 


Gi 


= X1X3X5X7 Gi 


= Z\ZzZ$Zi 


G 2 


= X^X^XfiXi G-2 


= Z^Z^Z^Zi 


G z 


— X2X3XSX7 Gs 


= ZiZzZ^Zi 



and so corrected. As seen, every multi-qubit error, is, up 
to stabilizers, equivalent either to a singlc-qubit error or 
to one of three two-qubit errors, 
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